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A  REMARK  ON  RADIATIVE  DECAY  RATES  FOR  MOLECULES  NEAR  A  DIELECTRIC  SPHERE 


P.  T.  Leung,  Young  Sik  Kim*  and  Thomas  F.  George 
Departments  of  Chemistry  and  Physics  &  Astronomy 
239  Fronczak  Hall 

State  University  of  New  York  at  Buffalo 
Buffalo,  New  York  14260 


ABSTRACT 

A  model  for  radiative  decay  rates  of  molecules  near  a  dielectric  sphere  by 
Gersten  and  Nitzan  is  analyzed  and  compared  to  the  results  from  an  exact 
dynamical  theory. 


*  Also  affiliated  with:  Department  of  Chemistry 

University  of  Rochester 
Rochester,  New  York  14627 


In  a  series  of  papers,  we  have  carried  out  critical  assessments  of  the 


validity  of  the  static  image  theory  (IT)  as  applied  to  the  problem  of  the 


classical  decay  rates  of  molecules  near  a  dielectric  surface,^  where  the  role 


of  the  dynamical  energy  transfer  theory  (ET)  has  been  emphasized  for  the  cases 


12  3 

where  the  surface  is  flat,  spherical  and  of  shallow  roughness.  In  particular, 


we  have  shown  that  IT  can  be  very  inaccurate  for  highly-conducting  substrates  and 


that,  strictly  speaking,  the  radiative  rates  cannot  be  described  within  IT. 


However,  it  is  interesting  to  note  chat  in  a  series  of  paper,  Gersten  and  Nitzan 


(GN)  succeeded  in  incorporating  the  radiative  rates  into  IT  for  a  spheroidal 


substrate  via  a  simple,  approximate  model.'4  It  is  the  purpose  of  this 


communication  to  assess  the  accuracy  of  the  GN  model  for  the  radiative  rates, 


where  we  shall  consider  a  dielectric  sphere  as  the  substrate  throughout. 


The  total  decay  rate  can  be  expressed  as  the  sum  of  the  radiative  and 


nonradiative  transfer  rates, 


R  NR 
7-7  +1 
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where  S  -  —  Re(ExH*)  is  the  Povnting  vector,  a  is  the  conductivity  of  the 

OTT 

substrate,  W  is  the  average  energy  of  the  dipole,  dfl  is  the  differential  element 
of  the  solid  angle,  and  dr  is  the  differential  volume  element.  The  integral  in 
Eq.  (2)  is  over  a  surface  at  infinity  while  that  in  Eq.  (3)  is  over  the  volume  of 
the  dissipative  substrate.  For  a  dipole  above  a  sphere  of  dielectric  constant  t, 

the  solutions  for  the  exact  dynamical  fields  and  hence  the  decay  rates  have  been 
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reported  in  the  literature,  ’  where  previous  comparison  with  the  static  theory 

(IT)  shows  that  as  long  as  the  distance  (d)  between  the  molecule  and  the 

spherical  surface  is  much  smaller  than  the  emission  wavelength  (A)  and  that  the 

conductivity  of  the  sphere  is  not  too  high,  IT  should  give  results  for  the  fields 

very  close  to  those  from  ET .  ^  Hence,  we  expect  that  7^,  which  is  important  only 

for  close  distances  d,  can  be  treated  with  reasonable  accuracy  in  the  non- 

retarded  theory  as  in  GN's  approach.^  However,  it  is  not  difficult  to  see  that 

R 

if  one  substitutes  directly  the  fields  from  IT  into  Eq.  (2),  one  gets  7  0 

2—*  A  2 

since  lim  r  S*n  =  0  for  fields  from  IT.  Thus  it  appears  that  a  straightforward 

£*■ -HD 

application  of  IT  cannot  incorporate  the  radiative  transfer  for  the  decay 

£ 

problem.  Nevertheless,  GN  has  considered  the  following  model  to  calculate  7  . 

4 

By  calculating  the  total  induced  dipole  moment  for  the  sphere  from  IT 

:in-blvd'  •  <4> 

the  dipole-sphere  system  is  then  regarded  as  one  free  dipole  system  with  total 
moment  oscillating  with  the  same  frequency  as  the  original  molecular 

dipole  .  Hence,  because  of  this  "single-system  approximation",  the  possible 
interference  between  the  fields  in  the  radiation  zone  emitted  by  the  two  dipoles 
and  4^^  (which  are  actually  at  different  locations)  cannot  be  included  in  this 


i 
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approach.  The  radiative  rates  (normalized  with  respect  to  the  free -molecule  rate 
7n)  thus  obtained  are  then  given  by^ 


R 
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where  a  is  the  radius  of  the  sphere,  and  and  ]|  denote  the  radial  and  tangential 

orientation  of  the  dipole,  respectively.  On  the  other  hand,  the  results  from  the 

2  5 

exact  dynamical  theory  (ET)  can  be  expressed  as  ’ 


71  3  Y  ,Jn(y}  +  BnhnL  (y) 
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where  j  and  lv  are  the  spherical  Bessel  functions,  y  -  2x(a+d)/A,  and  the  Mie 

coefficients  A  and  B  may  be  found  in  Ref.  5. 
n  n  y 

We  have  performed  some  numerical  comparisons  of  the  results  in  Eqs.  (5)  and 
(6)  with  those  in  Eqs.  (7)  and  (8),  where  we  have  generally  restricted  ourselves 
to  cases  with  d  «  A  and  a  «  A,  for  which  both  the  IT  and  dipole  approximations 
should  be  good.  Moreover,  these  conditions  are  often  realizable  in  actual 
situations.  However,  this  still  leaves  no  restriction  on  the  comparison  between 
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d  and  a.  We  have  performed  calculations  for  a  Ag  sphere  and  for  the  emission 

frequency  set  at  3  eV  (i.e..  A  -  4133  A) ?  Figure  1  shows  the  results  from  the 

GN  model  (dashed  lines)  compared  with  the  exact  theory  (solid  lines)  for  a  -  400 

A,  from  which  we  observe  already  appreciable  discrepancies  for  the  GN  model  with 

such  a  realistic  size  of  the  sphere.  We  have  also  repeated  the  same  calculations 

for  a  =  100  A  where  the  GN  model  has  indeed  been  found  to  lie  close  to  the  exact 

theory.  Figure  2  shows  the  comparison  at  a  fixed  d  —  100  A  with  various  values 

of  a,  where  it  is  clearly  shown  that  the  GN  model  is  accurate  only  when  a/d  <  2.5 

in  this  situation.  The  discrepancy  arises  because  of  both  the  dipole  and  the 

"single-system"  approximation  made  in  the  GN  model.  This  latter  approximation 

leads  to  the  fact  that  the  oscillatory  behavior  of  7  /y^  versus  a/d  totally 

disappears  in  the  GN  model.  It  may  be  of  interest  to  mention  that  similar 

oscillating  behavior  has  been  observed  for  the  flat-surface  case  due  to  the 

interference  between  the  reflected  field  at  the  dipole  site  with  the  dipole 

8 

oscillator  itself.  Moreover,  it  has  also  been  acknowledged  by  GN  that  their 
results  cannot  yield  the  correct  asymptotic  results  for  a  flat  surface  in  the 
limit  a  -*  °°  with  a  fixed  value  of  d,^  whereas  it  has  been  noted  that  Eqs.  (7)  and 
(8)  can  indeed  lead  back  to  the  correct  flat- surface  limit.  We  have  further 
shown  that  these  discrepancies  do  not  arise  from  the  dipole  approximation  alone. 
In  Fig.  3  we  show  similar  calculations  as  in  Fig.  2  for  the  range  0  <  a/d  <  4 

with  d  fixed  at  50  A,  which  implies  that  a/A  <  0.048.  Discrepancies  are  obvious 
in  spite  of  the  fact  that  the  dipole  approximation  should  hold  for  such  a  small 
ratio  of  a/A.  Hence,  we  conclude  that  one  must  be  very  careful  about  the 
limitations  of  the  GN  model  when  it  is  applied  to  the  dipole -sphere  system.  In 
addition  to  the  conditions  a  «  A  and  d  «  A  as  noted  before,  the  value  of  a 
cannot  be  too  large  for  a  fixed  d.  Since  some  of  the  previous  calculations  do 
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not  exhibit  an  awareness  of  this  condition,  it  would  be  of  interest  to  re¬ 
examine  all  the  problems  on  surface -enhanced  Raman  scattering  and  other 
photochemical  processes  as  previously  treated  by  GN  adopting  their  model  for  che 
molecular  decay  rates.4 
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Figure  Captions 

1.  Comparison  of  the  GN  model  for  the  radiative  rates  with  the  exact  dynamical 
theory.  The  sphere  has  a  fixed  radius  with  a  -  400  A,  and  all  other 
information  is  given  in  the  text. 

2.  Same  as  Fig.  1,  except  that  a  is  varied  and  d  is  fixed  at  100  A.  Note  that 
the  scale  of  the  y-axis  is  not  uniform. 

3.  Same  as  Fig.  2,  except  that  d  is  fixed  at  50  A. 
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